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SUMMARY 


An analysis is given of the oscillating motion of a ballistic missile 
which upon entering the atmosphere is angularly misaligned with respect to 
the flight path. Idle history of the motion for some example missiles is 
discussed from the point of view of the effect of the motion on the aero- 
dynamic heating and loading. The miss distance at the target due to mis- 
alignment and to small accidental trim angles \is treated. The stability 
problem is also discussed for the case where the missile is tximbling prior 
to atmospheric entry. 


INTRODUCTION 


It is characteristic of long-range rockets that, because of the low 
efficiency of the propulsion system, the weight at take-off is large com- 
pared to the final weight after fuel is expended. Typically > a saving of 
1 pound in final weight can save of the order of 20 poxonds of initial 
weight and, as a result, strict attention must be ~given"an“'the design of 
rockets to keep design safety factors to a minimum. Thus the magnitude 
of the factors which principally influence the final weight must be known 
with as great acciiraey as possible. 

Two such factors are the aerodynamic loal experienced by the warhead 
unit as it descends through the atmosphere, which affects required struc- 
tural weight, and the aerodynamic heating experienced in the descent, which 
affects required coolant weight. Problems relating to the loading and 
heating of missiles during atmospheric entry have, of course^ been given 
considerable attention, both from a general point of view (e.g., ref. 1) 
and in detail for specific designs. In the usual treatment of the problem, 
however, the rather idealized case has been treated wherein the vehicle 
enters the atmosphere unyawed or unpitched with respect to the flight path 
and without angular velocity. If the vehicle enters the atmosphere in a 
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yawed or pitched ajbtitude, it will, during its oscillatory approach to the 
earth, he subjected to lateral forces in addition to the longitudinal 
forces due to deceleration. Moreover, the distribution of aerodynamic 
heating over the surface for the oscillating vehicle will differ from 
that for the vehicle if aligned with the flight path. Thus a question 
arises as to what extent the structural weight and the weight of coolant 
might be altered by the fact that the rocket upon entering the atmosphere 
is angularly misaligned and has angular velocity. 

The analysis of reference 2 provides an excellent basis for such a 
study. However, the results of that analysis are in a form which is not 
convenient for demonstrating the relative importance of the several fac- 
tors which are of principal interest to the loading and heating problems. 
It is the purpose of this paper to re-examine the motion analysis of 
reference 2 using some simplifications which were employed in reference 1 
in order to Indicate more clearly the salient features of the motion 
problem and its effect, in turn, on loading and heating.^ 


BASIC ASSUMPTIONS 


In the analysis to follow it will be assumed that the missile warhead 
which enters the atmosphere is rotationally symmetric so the misalignment 
angle may be considered as yaw or pitch or any vector combination thereof, 
and that the fineness ratio is sufficiently low and the Ifech nxmiber is 
sufficiently high that the pressure distribution is independent of Jfe.ch 
number for the Mach number range of importance (see ref. 5) • Thus, the 
rates of change of the respective aerodynamic force and moment coeffi- 
cients with a, d, and q are considered to be constants. The basic 
assimiption of the analysis of reference 2 is retained in the analysis to 
follow; namely, that the angular oscillations are small so that the sine 
of the angle of oscillation is the angle of oscillation in radians and 
the cosine is unity, and the drag coefficient is sensibly the same as it 
would be for the nonoscillating missile. In addition, the assimiptions 
employed in reference 1 are also employed herein; namely, that the accel- 
eration of gravity is constant with altitude, the flight path through the 
atmosphere is essentially a straight line, and the variation of air density 
is the exponential function 


^At the time of the writing of this report, it was discovered that 
an analysis (ref . 3) very nearly identical to the analysis given herein 
was made by Dr. George Solomon of the Ramo-Wooldridge Corporation which, 
imf ortunately, has been given only a very limited distribution. However, 
the reference 3 analysis ignores the important effect of plunging on the 
damping of the motion. This effect is included in this paper. It shotild 
also be noted that an analysis has been made in reference k, which although 
principally aimed at study of other features of the stability problem, 
employs a basic approach similar to that of the present report. 
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( 1 ) 


vherein and p are constants and y is the altitude measured from 
sea level. 

A complete list of symbols is given in Appendix A. 


ANALYSIS 


If the angular displacements are small, then the differential equation 
of angular motion with time as the independent variable may be written 


d^g 

dt2 


+ fl(t) 


dg 

dt 


+ f2(t)g 


0 


( 2 ) 



wherein the time -dependent coefficients are given by^ 


^This formulation is equivalent to equations (12) and (I 3 ) in ref- 
erence 2 except that in reference 2 the value of has been tacitly 
assumed to be zero, which is a justifiable assumption at high Mach num- 
bers. The quantity 0^^ is retained in equation (3a) to be consistent 
with NASA standard nomenclature (e.g., see ref. 6). 
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fi(t) = 


Pm 


(Cuiq^ + Cin^)pVAl^ 
21 


(3 a) 


f2(t) 


_d^ 

dt 



iHm 


C^pV^AZ 

21 


(3^>) 


where the angle of attack Is as indicated in the sketch. 

It is convenient, now, to rewrite equation (2) with altitude, y, 
rather than time, t, as the independent variable. To this end, it is 
noted that for the straight flight path assumed (0 = % ” constant) 


dt 


-V sin 0g 


(i^a) 


d^y 


df 


2 S ^ % = V sin20E 


dy dt 


dy 




Thus 


dot. _ do, dy 
dt “ dy dt 


dot, 


= - ^ = -V sin 0E ^ 


(5a) 


dt" 


a _ ^ ^ _ vSginSfl 


dra 


<iy 


V VA. V « ^ VJUVJU 

— sin‘=^0-r. — 


da 
® dy 


(5b) 


so that equation (2) becomes 


d^a 


da 


^ + f3(y) ^ + f.(y)a 


= 0 


( 6 ) 


wherein the altitude -dependent coefficients are 
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f 4 (y) 


t ( ) = (Gm^ + 

® ^ V 2m sin Qg ^ 21 sin 0g 

ei^A(dp/dy) Ci^Ap(dV/dy) Cn^Cj^p^k^l^ Cjn^pAZ 
2m sin 0g 2mV sin 0g iHm sin20g 21 sin20g 


(7a) 


(7b) 


It is shown in reference 1 that by use of the same/hasic assumptions 
as have been made in this report, the velocity may be expressed as 


V = 




ko -w 
■5“ e 


(8a) 


where 


ko 


CpPpA 
pm sin 0J; 


m 


Hence 


d-V/dy ^ p^ -py 
V 2 ® 


and from equation (1) 


8(c) 


dp 

ay 



(8d) 


Thus equations (7a) and ( 7 b) become, upon setting the square of the radius 
of gyration as = l/m. 


f3(y) = 


PQ^ 

2m sin 9 ^ 


% - C + 


(Oniq^ + 


(9a) 


f 4 (y) = 


PcA 

2m sin^0. 


Cjjj^Z 


E' 


+ Ci^P sin 0g 


g-py^. Po^A2 




■Oi^Od- 




2n 


.-2py 


(9b) 
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If we set 


Y = ^y 


then equation (6) may he written 


d^a 


dY‘ 


dY 


-Y da /, -Y , -zY\ 

2 + 2kie — + ( kge + kae Ja = 0 




wherein the constants are 


ki = 


Po^ 


4pm sin 0g 


PcA 






ks = 






4p2m2sin20„ L 


0)T 


In order to find a solution for the differential eq^uation 


kie 


-Y 


a = Tje 


so that 


-Y -Y 

da kre dt] -Y kie 

— = — - kxe e Ti 


dY 


dY 


and 


d^a 

dY^ 


= e 


kie ^ d^T) 
dY^ 


-Y kre"^ dTi / -Y k,e p -aY kie 

2k le e U 1 ^le e + ki e e 


Substitution in equation (lOb) then yields 


(10a) 

(10b) 

(10c) 

(1Gb) let 

( 11 ) 
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-Y 

+ kx)e + 


(ie - ki^)e 



= 0 


( 12 ) 


The rigoroiis solution of this equation is not necessary since the 

term 

(ks - ki^)e"^^ 


is, for the type of missiles to he considered later, small coii^ared to 

(ka + ki)e"^ 


particularly at the higher altitudes which are of principal concern. If 
the square tenn is omitted, the solution of equation ( 12 ) is known (see, 
e.g. , ref. T) 




“ C iJof 2 kg + kj 




2>/k2 + ki e 


(13) 


where Gi and Ca are constants of integration and, in accordance with 
Watson’s notation (ref. 8 ), the ftmctions Jq and Yq are the zero order 
Bessel functions of the first and second kind, respectively. 

Combining equations (lOa), (ll) , and (I3) gives finally 


a = 




+ ki e 



+ kx e 



( 14 ) 


If it is specified that on entering the atmosphere (y -> «>) the 
missile is misaligned by the angle oog but has no angular velocity, then 


Cx - Oig 
Ca = 0 


(15 a) 
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and equation (l 4 ) can be written 


g 

“e 



+ kx e 



(15b) 


In the cases of usual interest the quantity Isq is very much larger than 
kx, as will be shown later, so that one may use the approximation 


g _ kxe 




(15c) 


For large values of the argment, the approximation 



(15a) 


can be used, and the maximum or envelope value is thus 




(15e) 


In Appendix B are derived expressions for the angular velocity and 
angular acceleration as well as their maximum values during any cycle of 
oscillation. The approximate maxlmimi angular velocity per cycle is (see 
eq. (B5 c)) 



- M 

pVgSm 0^e 



(I 5 f) 


while the maximimi angular acceleration per cycle is 
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"d^(a/ag)" 
- dt2 . 


* 


3 

4 






E 




E' 


spy 

4 


(ki-ko)e 


-py 


(I5g) 


From equation (l 5 d) It can t>e foimd that one cycle of oscillation 
takes place during the altitude change 


2 it 

P ./ka e ^ 

The frequency of oscillation is thus 


- M 

$ ./ka e ^ V sin 

(*) = — — — 

2rt 




and using the velocity from equation (8a) one has 


(0 


pVg sin 0E 

2rt 



M 

2 


(16) 




The frequency is maximum at the altitude 


yi 


1 

P 


In ko 


(17a) 


which, as shown in reference 1, is the altitude at which the deceleration 
due to drag is maximum and the velocity is 



0,61 Vg 


(17b) 


The maximum frequency is 


^ax 



_ i 


'pVjje “^sin % 




(17c) 
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and the corresponding an^litude is from equation (15©) 



(I7d) 


For later use it is of importance to note that in the analysis of 
reference 1 it is shown that for turbulent flow the altitude for which 
the average heat -transfer rate is a maximum is 



(18a) 


when 


_ 1 

V2 - e‘ ® Vg = 0.72 Vg 


(I8b) 


The corresponding frequency is 


0)2 


1 

S / pV^e ^sin 

^ V iiT J 


(I8c) 


and the corresponding amplitude is 



(I8d) 


On the other hand, for laminar flow (which, it is expected, should 
he applicable at least for the stagnation point on the missile nose) , the 
altitude for maximum heating rate is then 


73 = ^ 2n(3ko) 


(19a) 
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when 


1 

V3 = e“ ® Vg = 0.85 Vg 


(19b) 


and 



Prom a loads point of view it is of value to note that from equa- 
tion (15g) the altitude for which the maximum angular acceleration 
occurs is 


74 = p 


I (1^0 


^1) 


(20a) 


and that the maximum acceleration is 


d^Ca/oE)' 


L dt® 







^ 3 

e - p^Vji^sin^0jj 


(20b) 


Also the normal force experienced is 



(21a) 


and substitution from equations (l), (8a), and (l5e) gives, for the maximum 
nonpal force per cycle 


N* 


%a°^EPoVE^^ Py^(ki- ko)e‘^^ 

2J JTs/^ 


(21b) 
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This force is a maximum at and has the value 


M, 


max 


W4* = 


2. >J It \l lat-s 


Il—i—T -II 

lU(ko-ki)J ^ J 


(21c) 


Finally, it is of interest to determine the order of magnitude of 
the drift due to og from the course the missile wotild have had if aligned 
with the flight path. If n is defined as the distance normal to the 
straight line trajectory (that the missile wo\ild have if og were zero), 
the lateral aeceleration is 


d^n 

dt2 



m Pm 


( 22 ) 


Use of the relations of equations (l), (8a), and (l5c) gives the lateral 
velocity at altitude y 


dn 

dt 


^Iq,”E^EPcA 
2m sin 



(23 a) 


With the substitution 


_H± 

= 2^/^ e 2 


then this may be written 


dn 

dt 


e 

^Lg°E^EPcA r 

4k23m sin 9g 


py 

2 



?iJo(Si)<i^i 


(23h) 


The lateral displacement at y = 0 is, by similar substitution. 
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n = 


2 n/^ Ao^gg 

^i<x,°epcA r 

2ka 8 in^ ©E 


c 4 


r' A^y 


CiJo(5i)4Si<iS ( 23 c) 


Defining 


F(ko,ki,k2) = 


2 





i 


4ka 


CiJo(Ci)<i5a'i^ (24a) 


then 


ClaO-EPeA , , , 

n - — ^ ' ' .Lj . F(ko,ki,k2) 


2p^ sin^Qg 


(24b) 


The integration of eq.uation (24a) may be performed in the special case for 
which kx = ko = 0, for then 


F(0,0,k2) 



£iJo(ti)i^i 


2 Jlss 





1 - Jq(2 \I ^z) 

— 



For values of ko and kx other them zero, an analytic solution of equa- 
tion (24a) is not known. However, Dr. William Mersman has determined 
values of the function by numerical Integration using an IBM 65O type 
digital computer. The computational procedure and results are given in 
Appendix C. Over the range of variables of interest in this paper, 

0 < ko < 10, -80 ^ kx < 0, 5x10'* < k2 < 8x10®, the integration Indicates 
that 

^(koAi^ie) - 1A2 

within 0.2 percent (see Appendix C, eq.. (G20)). 
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If the angle is considered as a yaw angle then the miss distance 

is the lateral or "deflection target error" 



2p^m sin^Qg 


F(ko,ki,k2) 


(25a) 


while if the angle Og is considered as a pitch error, then the miss 
distance is the longitudinal or "range target error" 



2 p^m sln®0g 


F(ko,ki,k2) 


(25b) 


Thus^ in the general case^ the area of miss is elliptical with the 
major axis in the range direction (except when the missile descent is 
vertical when the miss area is circular) . 


It is useful^ for comparison purposes^ to find the miss distance which 
restilts from an accidental trim angle ^ Orji. The differential equation ( 22 ) 
becomes in this case 


d^n 

dt^ 



Ci^g^V^pA 

2m 


(26 a) 


SO that 


dt “ 2m sin 0^ 4) ' 


or 


dn 

dt 


' ' ' I 1 e 
koPm sin \ 



and, in tiorn. 


n 


^I<x^PcA TYl 

koPm sin20g \ 



(26b) 
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which can he shown to give 
wherein 

Fi(ko) = ~ [^(¥) - 0 * 5 t 72 l 6 j ( 26 d) 


and lT(ko/2) is the exponential integral for which tabulated values are 
given in reference 9* Values of Fi(ko) have been computed for practical 
values of ko and are presented in table X. 

If the trim angle is a yaw angle then the miss distance is the 
deflection target error 




sin^0g 


Fi(ko) 


( 27 a) 


while if the trim angle is a pitch angle then the miss distance is the 
range target error 



- PlCko) 

sin®0-g 


DISCUSSION 


(2Tb) 


Conical Warheads 


It is the purpose in the discussion to follow to examine the angular 
motion of typical ballistic missile warheads in the atmosphere to ascertain 
the importance of this motion to the problems of aerodynamic loading and 
heating and miss distance. Conical shapes have been chosen for this study 
for the reasons that, first, they are often considered suitable for appli- 
cation to warheads and, second, the calculated stability derivatives are 
available . 

Tobak and Wehrend (ref. 10) have calculated the stability derivatives 
for cones of half- angle, 6. Althoxigh they give results which are applicable 
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from low supersonie speeds (exceeding the Mach number for shock detachment) 
to hypersonic speeds, the concern of this paper will be only with the 
hypersonic or Newtonian solutions. For eirbitrary distance from the cone 
apex to the center of gravity, leg, reference 10 gives the pertinent 
derivatives, using the symbols of this report, as 


Cw = 2 cos% 

Cx = Cj) = 2 sin^S 




= ^N„, - Cp = 2 (cos 26 - sin^B) 

= -(1 + tan^B) + I - 2 cos^b^Y^^ 








= 0 


- — + 2 cos^B 



J 


(28a) 


In addition it should be noted that for cones the center of volume is at 
3Z/^ and the square of the ratio of cone length to radius of gyration for 
arbitrary center-of -gravity position is (considering the body to be 
uniformly solid) 



80 


12 tan^6 +3+80 



(28b) 


Validity of the Analysis 


Before examining the stability, loading, heatirg, and miss-distance 
problems it is necessary to detemine whether the previous analysis is 
valid for the conical shapes to be considered. It was noted that the 
solution given by equation (15b) applies only in those cases for which in 
the differential equation (12) the term 
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(ks - (29a) . 

can be neglected in comparison with 



(29b) 


In Appendix D it is shown that the values of the factor (29a) are, 
for practical cases, always very small compared to the Values of the 
factor (29b) and, in addition, that ka is, for practical, cases, always 
much larger than kx so that the solution of equation (l5c) is valid. 


Stability 

In the equation (l5c) it is dear that the missile is statically 
stable if kg is positive and certainly dynamically stable if kx is zero 
or less than ^ero; If kx is positive, then it is possible for the oscil- 
lation amplitude to increase with decreasing altitude. This may be con- 
veniently shown from the approximate equation (l5e) for the envelope value 
of a. Writing this equation in the form 



we find that the derivative of this function with respect to y is zero 
when 


1 "3y 1 

kxe 


or 


y ^ I ln(4kx) 


(31) 


Thus if positive kx is even as large as 0,25 the incipience of divergence 
occurs (at sea level) . For larger positive values serious divergence at 
the lower altitudes must be anticipated. In order to obtain a better grasp 
of the nature of these motions it is instructive to examine the behavior 
of a missile with some arbitrarily assumed (but practicably realizable) 
static stability for several values of kx. For this purpose let it be 
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supposed tliat kg is 10'^ and that ki has, in turn, the value -10, 0, 
and +10. The equation (l5e) then yields the envelope curve of angular 
history shown in figure 1. j^t is clear that positive values of ki could 
promote serious divergence of the amplitude of angular oscillation. If 
a missile had a large positive value of k^, as shown, it would almost 
certainly tumble at the lower altitudes. Friedrich and Dore (ref. 2) dis- 
cussed the possihllity of such divergence and noted that these adverse 
effects could occur if the missile underwent large reductions in speed 
during the descent. ^However, no consideration was given in their report 
to the importance of the damping terms* In other words what they con- 
sidered, in the language of this report. Was that the coefficient ki (as 
given by eq. (lOc)) was overwhelmingly Influenced by the drag coefficient. 
■Hiese results of their analysis thus imply that high drag shapes are unsat- 
isfactory for ballistic missile application in spite of their inherent 
advantages in the aerodynamic heating problem (see, e.g., ref. l). 


The question of importance is, then, 
to ignore the damping terms (Cj^) and (Cm 


whether or not it is realistic 
+ Cm^) due to plunging and 


rotation, respectively, in the determination of sign and magnitude of 
ki. To answer this question it is convenient to consider a simple conical 
warhead of arbitrary cone half-angle, 8. 


To investigate the sign of k^ for cones it is sufficient to evaluate 
the "dsmamic stability" factor 


^D - + (Cmq + 

in equation (lOc) by use of the relations of equations (28a) and (28b). 
This has been done and the results are presented in figure 2 for several 
center-of -gravity positions. (Note that Iq^/i = 3 A ts the center of 
volume and a most likely position for the center of gravity. ) From fig- 
ure 2 it is seen that ki must always be negative for conical (and pre- 
sumably for near conical) shapes. Thus the inference of reference 2 that 
high drag (i.e. , large 8) shapes are undesirable from the dynamic 
stability viewpoint is unjustified. 

Next it is in order to exam i ne the sign of the factor which controls 
the static stability parameter 


It is shown in Appendix D that the second term in the bracketed expression 
is. small compared to the first for practical cases so that it is only 
necessary to be assured that is negative to Insure stability. In 
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figure 3 "the Newtonian value of (from eq.s. (28a) and (28b)) is 

plotted as a fimction of 6 for various values of Icg/l, and it is seen 
that when 6 is small, care must he exercised to keep the center of gravity 
far forward. It shoiild he noted (see ref. 10 ) that the center of pressure 
is independent of tfe,ch number down to the Mach number of shock detachment. 
Thus the hypersonic requirement of Z^g is also the supersonic requirement. 

From the preceding discussion, it is apparent that at least from low 
supersonic to hypersonic speeds, positive static stability (i.e., positive 
kg) can be obtained. Similarly, it can be shown that over the same speed 
range dynamic stability is assured (i.e., negative kj.) . Now it is impor- 
tant to determine the magnitude of the static and dynamic stability which 
can be provided. To this end the following digression in the discussion 
is in order. 


For long-range ballistic missiles the aerodynamic heating problem 
must be the principal consideration in design. It has been shown (see 
refs. 1 and 11) that, generally, the aerodynamic heating problems are 
reduced when the value of ko is increased. On the other hand, if ko 
is too large then the speed of descent of the missile becomes low for too 
great a part of the final trajectory which increases the vulnerability of 
the missile. Thus some compromise is required and this compromise value 
of ko tends to be larger the longer the range. A value that will be 
considered herein to be a reasonable one for a 3^000- to range 

would be of the order of 5 bo 20. 

Now, kx (see eqs. (8b) and (lOc)) can be written in terms of ko in 
the form 



while for kg, if the sin part is neglected (see Appendix D), this 

parameter becomes 


ka 


2pZ sin % 



(32b) 


Values of kx/ko and PZ sin ©gka/ko are given in tables II and III. Since 
kx/ko depends upon the location of the center of gravity and the cone 
angle while ka/ko depends upon these factors and, in addition, the length 
of the missile, it is necessary to consider some examples in order to 
determine the magnitudes likely to be realistic for kx and ka. Accord- 
ingly, let it be assumed that the missile weight is 3 >000 pounds, that the 
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entrance angle of the trajectory, Is 30°^ that the values of the 
atmospheric density relations sire those of reference 1 ( pq ,= 0.003^ slugs 
per cubic foot, 3”^ = 22,000 feet). For values of ko equal., in turn, 
to 5, 10, and 20 the base diameters, ler^ths, and volumes of the example 
missiles are those given in figures ^(a), ^(b) , and ^(e), respectively, 
as functions of cone half-angle. For this analysis it is arbitrarily 
assumed that the maximm allowable missile length is 30 feet and the mini- 
mum allowable volume is 10 cubic feet (corresponding to a high missile 
density of 3OO poionds per cubic foot) . Ihus in figures k the curves extend 
only to the cone half- angles which correspond to these two limits (the 
small cone-angle limit corresponds to the maximum allowed length and the 
large cone- angle limit to the maximum allowed density) . In addition, it 
is arbitrarily specified that the center of gravity in each case is at a 
distance from the apex (Zcg) where the local diameter is 2-1/2 feet. The 
resulting ratios of Zcg/^ shown in figure ^(d) , 

With these physical characteristics, the values of kx and kg are 
those of figures 5 and 6. It is seen that the dynamic stability is great- 
est for large values of the drag parameter but for small values of the cone 
angle. On the other hand, the static stability is generally greatest for 
large values of both the drag paraineter and cone angle. A notable excep- 
tion to this trend of the static stability parameter is the sharp decrease 
of k2 at the largest angle for the ko *= 5 case. This sudden reduction 
results from the fact that the center-of-gravity position has rather 
closely approached the center of press\ire. 

Consider, now, two extreme cases: first, the ko - 5 missile at maxi- 

mum allowed density and second, the ko = 20 missile at maximum allowed 
length. The former has least values for both stability parameters and 
therefore will oscillate with the largest amplitudes, while the latter has 
the largest dynamic stability parameter and has a rather high value for 
the static stability parameter and thus shoilLd be representative of the 
opposite extreme. The angular behavior with altitude for these two mis- 
siles has been calculated using equations (l5c) and (l5e) and is shown in 
figures (7a) and (7b). The high altitude oscillations of fig^lre 7(a) are 
similar to those of (7b) but displaced downward, altitudewise , by about 
25,000 feet. This is an effect of the lower static stability parameter 
for the ko = 5 missile. At the lower altitudes the ko = 20 missile 
oscillations decrease more rapidly by virtue of the larger dynamic stability 
parameter. 


Heating 


It was pointed out in the Introduction of this report that when the 
time rates of eierodynamic heating are largest, it is important that the 
oscillation amplitudes be small in order that additional coolant mass will 
not be required to protect the vehicle from excessive local heating. The 
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experimental data of reference 12 indicate that, at least for not-too- 
small cone angles, maximum oscillation amplitudes of the order of a few 
degrees of are shoiild he permissihle with no important adverse effects. 

To determine whether or not the oscillations will he in^jortant as 
regards aerodynamic heating it is again convenient to consider a particular 
example. The same 3 , 000-pound missiles are used in this study. Since the 
laminar heating rate always reaches a maximum at a higher altitude than 
does the maximum for turhtilent heating, it follows that (a/a^)* at maximum 
heating will he greater in the laminar case. The amplitude ratio at the 
altitude f or maxim\mi laminar heating rate (calculated using eg.. (I 9 d) ) is 
shown for the example missiles in figure 8 . It is seen that these values 
are so low that no complications of the maximxm heating rate prohlem due 
to Initial angular misaligranent with the flight path, otg, should he 
expected. In figure 9 hoth the angular amplitude ratio and the ratio of 
laminar heating rate to maximum laminar heating rate (see ref. 1 ) are 
plotted as a fiinction of altitude for one particular example (ko = 5 # 

6 = 25 °), which shows that while the aagjlitude ratio is very small at the 
altitude for which maximum heating occurs , it may hecome sufficiently 
important at the higher altitudes where the heating rate is still fairly 
high to require consideration in design. 


Loads 


To show the degree to which latereQ. loads due to 0 ^ are important 
it is again usefiil to consider the example missiles considered estrlier. 

To evaluate the maximum nomnal force using equation (21c) it is necessary 
to specify the entrance speed, and it will he assumed, for the exam- 
ples, that the speed is 20,000 feet per second. The maximum normal forces 
in terms of missile weight per degree angle misalignment at atmospheric 
entrance for the examples are shown in figure 10. While the normal forces 
are increased for the longer missiles (due to increased surface area) , it 
does not appeeir that they could be too serious in a practical case. A 
20 ° value for aj. only promotes a 3 g normal acceleration for the long 
missiles which is small compared to the deceleration due to drag which 
(from the analysis of ref. 1) is 51g* Moreover, the maximum nonnal loads 
are not additive to the maximum drag loads since, as seen in figure 11 , 
they occur at different altitudes. 


Miss Distance 


Before considering the actual magnitudes of the miss distances due ■ 
to Ojj or dip it is well to discuss the accuracy of the analysis of miss 
distance given previously. In the analysis it is assumed that the velocity 
at all points of the trajectory is given by the exponential ejcpression of 
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equation (8a) and tMs expression ms, in turn, obtained "by neglecting the 
effect of gravity. Hie neglect of the effect of gravity is unimportant in 
the evaluation of the velocity-altitude history except when the velocity 
is low and, simultaneously, the deceleration due to drag becomes comparable 
to the acceleration of gravity. Then the velocity given by the analysis 
falls below that which would actually occur. In the analysis of miss 
distance it should be clear that the miss distance increases rapidly as 
the velocity decreases, and, hence, the miss distances given by the 
analysis are in error by an amoiint which increases rapidly with increasing 
ko when the decelerations near y = 0 become of the order of the accel- 
eration of gravity. To assure that the deceleration at sea level is not 
less than Ig, it is required (from ref. l) that 


koe 


-ko < 


For example, for an entrance speed of 20,000 feet per second and for 
= 30°, the sea level deceleration reaches Ig for ko of about 7 (see 
ref. 1)5 hence, the values of F(ko,k2.,k2) (see Appendix O) and Fx(ko) 
from table I shoiild not be used, under these conditions, for values of 
ko in excess of about "J , partietilarly as ko greatly exceeds this value. 
Thus in the calciilations to follow the miss distances for ko = 20 are not 
included and, in addition, the reader must note that even for ko = 10 the 
estimated miss distances exceed the actual ones. 

The range target errors per degree angle misalignment at atmospheric 
entrance for the example missiles (Vg = 20,000 ft/sec, 0-g = 30°) are shown 
in figure 12 . The deflection target error, not shown, is simply one-half 
the range value. It is seen that, as with the normal force, the miss 
distance is greatest for the smallest cone angles. However, the miss dis- 
tance is trivial since even for an Og of say 20° the range target error 
is but about 20 feet in the worst case. 

A serious problem is the miss distance which will result from a trim 
angle even slightly different than zero. In figure I3 is shown the range 
target error per degree of trim angle for the ko of 5 and 10 . It is 
seen that a trim angle of as little as 0 . 1 ° can cause a range error of 
many miles. As the cone angle increases then for a given value of ko, 
the miss distance diminishes ■until when the cone half-angle is ^5° the 
miss distance is zero since the lift-curve slope is then zero (see 
eq. (27a)). Except in this special case, however, the miss distance due 
to even a slight trim angle is very important. One method for reducing 
this miss distance would be to spin the missile about its axis so that 
it would follow a corkscrew path diuring descent. This solution introduces 
another difficulty, however, in that care may have to be taken to keep the 
spin rate from approaching the pitching (or yawing) frequency else tiambling 
may occur if the missile is not identical as regards aerodynamic and 
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inertial characteristics about any raxiial axis (see ref. 6 or I3). 
Unfortunately, the pitching frequency varies from zero to the maximum 
value given by equation (I6) . The maximum frequeixcies for the example mis- 
siles are shown in figure l 4 . One obvious way to avoid tumbling resulting 
from "roll coupling" would be to spin the missile at a rate which exceeds, 
by a good margin, the maximum shown in figure 1 ^ 1 -, If the spin rates were 
the maximum pitch rates , the rim speeds at the base (i.e . , at maximum 
diameter) would be those shown in figure I5. Since the required spin rate 
would have to materially exceed this rate, it is clear that a serious 
stress problem due to centrifugal loading might result (especially for 
small cone angles). In consequence, a better solution might be to spin 
the missile at a rate which is always less than the value of the pitch 
frequency at any given altitude, but the difficulty, then, would be one , 
of assuring that the spin rate could not accidentally approach the pitch 
frequency. In the preceding discussion it has been tacitly assumed that 
the missile is not identical as regards the aerodynamic and Inertial char- 
acteristics about any radial axis so that the spin rate must not at any 
time match the pitch rate. Since the pitch rate changes rapidly with time 
(particularly at the higher altitudes), it is probably not a Justifiable 
requirement that the spin rate not ever be the pitch rate, particularly 
since the asymmetries which exist may be trivial. Some further eonsider- 
ation must clearly be given this problem. 


Effect of Initial Tumbling 


In the discussion to this point it has been assumed that the missile 
enters the atmosphere misaligned by an arbitrary but fixed angle with 
respect to the flight path. When the missile is actually tumbling before 
entering the atmosphere, then the analysis given previously cannot be used 
since the equation of motion is restricted to small-angle considerations. 
(That the analysis is inadmissible is reflected in the fact that if, in 


equation (l 3 )^ C2 


has a value which is not zero, then 



beccmes infinite if y is infinite.) In spite of this deficiency, some 
general remarks can be made about the effect of initial t\mibling. It is 
clear at the outset that the missile must have but one possible trim 
attitude if initial tumbling occurs. If not, it could descend at some 
attitude for which no adequate protection for aerodynamic heating and 
loading had been provided. Furthermore it is a requirement that it must 
be righted to about the correct attitude at an altitude which is suffi- 
ciently high that the angular motions will become small by the time the 
heating and loading are intense. One obvious way which might assmre that 
these conditions will be met would be to leave the empty fuel and oxidizer 
tankage connected to the warhead until the missile has entered the atmos- 
phere sufficiently fax to adequately correct the attitude, and then release 
of the tank could be pezmaitted. 
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CONCLUSIOnS 


From an analysis of the motion of a halllstlc missile Initially mis- 
aligned with respect to the flight path prior to the entry Into the 
atmosphere. It Is concluded that It Is posslhle to; 

1. Provide a continuously damped oscillation history with descent 
through the atmosphere. 

2. Keep the oscillations to a small amplitude when at altitudes for 
which aerodynamic heating and loading are severe. 

3 . Prevent excessively large loads due to the oscillating motion. 

Moreover, while the miss distance at the target due to the Initial 
misalignment angle Is trivial, the error that can occur due to the trim 
angle being even slightly different than zero can be very large and Its 
effect must be minimized In some manner. 

Since tumbling may occur prior to entry Into the atmosphere the mis- 
sile must have only one trim attitude and must be brought near this attitude 
before the missile has progressed too far down through the atmosphere. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., June I 5 , 1956 
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APPENDIX A 
SYMBOLS 


constant (See Appendix C, eqs. (C2),) 

reference area for coefficient evaluation (base area 
for cones) 

constant (See Appendix C, eqs. (C2).) 
constant (See Appendix C, eq.s, (C2).) 
constants of integration 
drag coefficient 

rate of change of lift' coefficient with angle of 

attack , f ) 

\da /a ^ 0 

rate of change of mcment coefficient with angle of 
attack . 




a -> 0 

change of moment coefficient with time rate of change 
of angle of attack. 


ft 


d 0 

rate of change of moment coefficient with angular 
velocity. 


(• 0 . 


q 0 

rate of change of noimial- force coefficient with angle 
SCn 

of attack, — — — 

0 

axial- force coefficient 
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d 

e 

fsCy) ^f4(y) 

Fi(ko) 


P(ko,ki,k2) 


g 

h 

I 


^0,1,2,.. . 

Jo( ) 

Ji,2,3,.*.,r,s( ) 
ko 


ki 


k2 

ks 


2 



m 

n 
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diameter of body base 
Naperian base 
fxmctions of time 
functions of altitude 

function used in evaluating of "miss distance" due 
to o<ji (See table I.) . 

function used in evaluation of "miss distance" due 
to (See Appendix C.) 

acceleratipn due to gravity 

integer (See Appendix C.) 

mass moment of inertia 

integrals (See Appendix C. ) 

Bessel function of the first kind of zero order 


Bessel function of the first kind of order 1,2^3^ • • • ^3 t,s 

the "drag" parameter (See eq.. (8b).) 

the "dynamic stability" parameter (See eq. (lOc),) 

the "static stability" parameter (See eq. (lOc).) 

the "cross-products" parameter (See eq, (lOc).) 

body length and reference length fol* moment coefficient 
evaluation 

distance from body bow to center of gravity 
cross-wind force 
missile mass 

distance normal to the trajectory the missile would 
have if it were angularly aligned with flight path 
and without angiilar velocity 

normal force (force perpendicular to the axis of 
revolution) 
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functions (See Appendix G.) 

<1 

angular velocity 

r 

Integer (See Appendix C.) 

s 

Integer (See Appendix C . ) 


functions (See Appendix C.) 

t 

time 

V 

speed at arbitrary altitude 

Vg 

speed on entry to the atmosphere 

X 

along range distance 

y 

altitude 

y 

dimensionless altitude, py 

yo( ) 

Bessel function of second kind of zero order 

z 

across range distance 

a 

angle of attack 

a® 

angle of attack on entry to the atmosphere 


angle of trim 

3 

density exponential (See eq. (l),) 

6 

half- angle of cone 

Kl 

deflection target error due to 


range target error due to og 


deflection target error due to 


range target error due to aip 

i 

altitude variable (See eq. (23b).) 


angle -of- attack function (See eq. (ll).) 
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P 

pQ 

a 

w 


angle between fllglit path and earth's surface that 
missile has on entering the atmosphere 

air density 

air density at sea level 
radius of gyration 

an arhitrary variahle (See Appendix B.) 
oscillation frequency 


( )* 


maximum value of the bracketed parameter which occurs 
in any particular cycle of oscillation 


Suhscripts 

( ) maximum value of the bracketed parameter 

' 'max 

Except for the parameters C, f, F, and k: 

( ) ^ value of the bracketed parameter at altitude for 

maximum deceleration 

( )g value of the bracketed parameter at altitude for 

maximum turbulent heat- transfer rate 

( )g value of the bracketed parameter at altitude for 

maximxjm laminar he at- transfer rate 

( value of the bracketed parameter at altitude for 

maximum normal force 
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APPEllDIX B 

DETEEMIMTION OP AEGUIAE VELOCITY AND AKGUIAE ACCEIEEATION 


It is the piirpose in this appendix to derive the ejcpressions for the 
angular velocity and acceleration from equation (15b) which gives the 
angular displacement as a function of altitude. Noting that 


dJo(t) 

d-y 


-J ) 


dt|f 

<iy 


then differentiation of equation (15h) yields 




ay 


1 

2 


•J^2 ^/k2 + e ^ + 
r _ M / Py\ 

\2s[kr+%. e 2 e J\ (Bi) 


and further noting that 


UJi(t)3 = tJo(t) 


dT|f 

ay 


then 


= ^2e^^® ^^(ki2e"2Py - kae"^^) jJ^ e" ^ - 

s3y , -py / J;A 

2p2kl^/k2 + ki e ^ Ji(2 e 2 J 


(B 2 ) 


Using equations ( 5 a), (8a), ajad (Bl) 
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APEEMDIX C 

NUMEEICAL INTEGEATION METHOD FOR THE lETEGRAL FCko^ki^ka) 


The following method for the numerical solution of the function 
F(ko,k i^ka) was devised by Dr. William Msrsman of Ames Aeronautical 
Laboratory. The integral to be evaluated is , from equation (24a) , 


pa 

kaFCkoAx.ka) = J ^ Io(^,c)d^ 


(Cl) 


where 


and 


a = 2-Jl^ 
b = ko/8k2 


c - Itt ko *" 5s:iy/4k2 




io(t,c) = 


(C2) 


(C3) 


By reference 8, page 

E ”’'Jr(Ci) = 3^ [?i”\+i(Ei)] . I- = °A,2,. . . (cl.) 


Introducing 


Ir(C,c) = , r = 0,1,2,. . . (C5) 


an integration by parts gives the recursion formula 
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That is, 

Ir(?,c) = + 2clr+i(^,c) , r = 0,1,2 (c6) 


/ \ X*"" s 

If this equation is multiplied "by (2c) and then summed on r 
from s to infinity, the following series representation is obtained 


Is(?,=) = y (2o£)>'J3^h^^(?) (C7) 

h=o 


for any s = 0,1,2,. . ., the series being convergent provided that 

|2cC| < |. 

In particular, then. Setting s =? 0 gives 


Io(?rc) = Y, (2cU^Jh+i(0 

h=o 


and substitution in equation (Cl) gives 


k2F(ko,ki,k2) - Y Ss(a,b,c) 
s=o 


(C8) 


where 


Ss(a,b,c) = (2c)® J" , s = 0,1,2, , 


(C9) 


the series is convergent if | 2ca[ < | . 
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Consider, first, the term s = 0 


So(a,b,c) = J' 


By reference 8, page I 8 , Ji(^) = -dJo(^)/d^. Hence integration hy parts 
gives 


So(a,h,c) = Jo(0) - e"(°'^^®-^Jo(a) - 2(c-h) 


Referring to equation (C3 ) , this is, since Jo(0) = 1, 

So(a,h,c) = 1 - - 2(c-h)lo(a,c-h) 


(CIO) 


To obtain a recursion formula for the general Sg, substitute in 
equation (C9) the equation 


Js+1«) = f Js(S) - 


from reference 8, page 45. This gives immediately 

Sg(a,b,c) = 4csSs_3^(a,b,c) - (2c)®Ig_3_(a,c-b) , s = 1,2,3,. . . (Cll) 


Thus, equations (c8), (CIO), and (Cll) reduce the problem to one of com- 
puting the sequence Ig(a,c-b), s = 0,1,2,. . . . 

For computing purposes, it is desirable to introduce slowly varying 
quantities. The following substitution turns out to be convenient 

Pr(a,c-b) = e'^^'^"^^®'\p(a,e-b)/a^+^ r = 0,1,2,. . . (C12) 


The computing problem is then summarized by the following formulae; 
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Pg(a,c-b) 


^ [2a(c-b)]^Js^h+x(a) 
li=o 


(C13) 


for any s = 0,1,2,. , . . 


Pi.(a,c-b) = Jr+i(a) + 2a(c-‘b)Pp+i(a,c-'b) , r - 0,1,2,. . . (Clil-) 
So(a,b,c) = 1 - e [Jo(a) + 2a(c-b)Po(a,c-b) ] (CI5) 

Ss(a,b,c) = l^csSs_i(a,b,e) - (2ca)V^°“^^®’Ss.i(a,c-b) , s =1,2,3,. . 


k2P(ko,ki,k2) 


00 

z 

s=o 


Ss(a,b,c) 


(CI6) 

(C8) 


The order in which the computations are to be performed is dictated 
by the following inequalities, each of which is obvious from the corre- 
spondixig integral definition, under the following general assumptions 


b > 0 , c > 0 , c-b > 0 
2ca < 1 , 2a(c-b) < 1 


The inequalities are 

lPs(a,c-b)| <— — , s = 0,1,2,. . . 

1 - 2a(c-b) 

|Ss(a,b,c)| <— (2ac)® , s = 0,1,2 
' * S -h 1 


(C17) 


(018) 


From the latter, it can be determined how many terms of the series (C8) 
are needed for any desired acc\xracy, so that (C8) is replaced in 
practice by 
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K 

k2F(korki,k2) = ^ Ss(a,-b,c) 

s==o 


(C19) 


In the present work K = 20 , giving a truncation error of less than 10 “^°. 
Once K is chosen, Pg- is obtained from equation (C 13 ) with s = K where 
again enough teims are taken to instire the desired accuracy. In the pres- 
ent work the con5)uter autcanatically continued the series (C 13 ) until the 

stonmand [ 2 a(c-h) became less than 10 ”^. Once % has been 

obtained from the series (CI3) , the recwrence relations (Cl 4 ) are used 
to con^)ute Pg-i, ^K-sf • • that order. Then Sq is confuted 

from equation (CI5), and Sx, S2,. . ., Sg in that order from 
equation (Cl 6 ) . 

In the present paper the significant range of parameters is 


G ^ ko 10 
-80 < ki < 0 
5x10^ < k2 < 8x10® 


For this range it can be shown that the series (GI9) can be truncated at 
K = 1 with an error less than 1 x 10 "®. Furthermore, in the expression for 
So and Sx, equations (CI5) and (Cl 6 ), the terms involving the exponential 
function are also less than 1 x 10 "®. This gives the simple approximate 
formulae 


So - 1 , Sx - 


and, hence 


F(ko,kx,k2) 



I ko " >1' 
+ — 



(C 20 ) 


with an error less than 2xl0"®/k2. 
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APPEHDIX D 

ORDER OF MAGNITUDE OF FACTORS AFFECTING STABILITY 


In the analysis of this report, a nnmher of simplifying assnD^jtions 
were made regarding the relative importance of the several factors which 
influence the stahility. It is the p\irpose herein to demonstrate that 
three of the assumptions which are of particular importance are, in fact, 
justified. These assumptions are: 

(1) In the evaluation of values of kg (sq. (lOc)) 
it is permissihle to ignore the Cj^ contribution 

in comparison with the Cm^^ contribution, (This 
assumption is desired but not required.) 

(2) In terms involving ka + k^, that ki is 
unimportant. (This assumption is desired but not 
required . ) 

( 3) That ks - kj_^ is trivial in comparison with 
k2 + kx (or, from assumption ( 2 ), in comparison with 
kg) . (This assumption is required to obtain the 
solution for the fimdamental differential equation of 
motion (eq. (l 2 ),) 

It is to be noted at the Outset that kx, and kg can be written 
in terms of the drag parameter, kg, in the forms 



The demonstration of the validity of the assumptions (l), (2), and (3) 
will be considered in the sections I, II, and III as follow for conical 
shapes. 


I. For this demonstration it is necessary to show that 
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sin 0g 


is trivial in comparison with 


Gd W 


Since the term is larger the larger the value of I and sin 9-^j 

it will he assumed that the is about (slightly greater than) 

and the length is that for a 6 -foot diameter base (i.e., I = ( 3 /tan 5) ft). 

Then (since = 22,000 ft) the comparison of the components is between, 

approximately 


and 


ClgXJ-O'* 

Cptan 6 


n / \2 

A ) 

Gd W 


Assuming I eg/ 1 is O. 3 O, O. 5 O, and 0,70 then the ratio of the exact value 
of kg to the approximate value of kg obtained by ignoring the 

term is that shown in figure I 6 , It is seen that the approximation is 
excellent except when goes to zero (shown for the ^cg “ O. 7 O case). 

Of course, in no practical case would Cu]^ be allowed to approach zero. 


II. To show that ki 
necessary to show that 


is trivial in comparison wit,h kg it is 







is small compared to 

^ ^ . 1 r Af ’ 

ko “ ‘ 2pz sin Sj. [' C-q \oJ ^ 
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Again the quantity kg will he least cong)ared to ki when I sin Qg 
is largest; hence, the length and angle assumption of section I is used 
so that the cQii 5 )artson will he between 

and 

kg _ tan 5 j\ 

ko ' 2X10-'^ \ W W 

Assuming ^cg/^ = 0*30, O. 5 O, and O. 7 O then the ratio of the value 
kg + ki to the value kg is that shown in figure I 7 . It is seen again 
that, except when Cju^ goes to zero, the approximation is excellent, 

III. In this demonstration it is necessary to show that 


(ks - kl 2 )e- 2 ^y = ko= 


is trivial con^iared to 



-2py 


(kg + ki)e-Py S kge-Py 



The test of this assxm 5 )tion is more severe when ko is largest and y 
is minimum, and again, when I sin 0 ji is largest. Using the largest 
ko to be expected, say 20, and the length and angle from section I, and 
for y = 0 , the severe cong)arison is then between 
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Assumiijg 1 q^/i - O.3O, O.5O, and 0,70, then the ratio of (l^-kx^)/k 2 
is that shown in figure 18. It is seen again that, except when Cnijx 
goes to zero, the approximation is excellent. Moreover since these coef- 
ficients enter the differential equation (l 2 ) multiplied by e'^Py and 
e"Py, respectively, and these exponential values have, at altitude, the 
values 


y^ 

ft 

e-2gy 

e" 3 y 

0 

50.000 

100.000 

150.000 

200.000 

250.000 

1.00 

1.06X10“^ 

1.12X10"'* 

1.19X10"® 

1.25x10"® 

1.31x10"^® 

1.00 

1.03x10"^ 

1.0^10"^ 

1.09x10"® 

1.12X10"'* 

1.15x10"® 


o 

then the integrated influence of the term involving k3 - ki must he 
trivial compared to the kg term. 
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TABLE I.- ERROR FTMCTION, Fi(ko), EQUATION (26d) 


ko 

Pi(ko) 

0 

0.500 

1 

•571 

2 

• 659 

3 

.773 

k 

.921 

5 

1.116 

6 

1.376 

7 

1.728 

8 

2.208 

9 

2.873 

10 

- - --- ^ 

3.800 


Note: The values should not he used when kg is such that 


S 


2g 

Vg p sin 
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Figure [.-Effect of parameter k, on angular amplitude ratio. 
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Figure 2.- Dynamic stability factors for conical missiles. 
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Cone half -angle , 8 , deg 


(a) Body base diameter. 

Figure 4,- Physical characteristics of example missiles. 
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Cone half- angle, 8 , deg 
(b) Body length. 

Figure 4.- Continued. 


COHFIDEHTIAL 




50 


C01PT3)11STIAL 


MCA EM A56FI5 



Cone haif-angie , 8 , deg 
(c) Body volume. 

Figure 4.- Continued. 
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Dynamic stability parameter 



Cone half- angle, 8, deg 

Figure 5. "-Dynamic stability parameter for example conical 

missiles. 
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Figure 8. - Amplitude ratio at altitude for maximum laminar 

heating rate. 
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Figure 9.- Laminar heating rate and amplitude ratio as 
a function of altitude for ko-5 , 8 = 25°. 


CONFIDEISTIAL 
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Figure 12. - Range target error per degree initial angular 
misalignment for example missiles. 
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Cone half- angle, S , deg 

Figure 13.- Range target error per degree trim angle for 

example missiles. 
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Figure 14. - Maximum oscillatioii frequencies for example 

missiles. 
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Cone half- angle, 8 , deg 


Figure 15.- Rim speed for example missiles If spin rate 
equals maximum oscillation frequency. 
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Figure IT. -Ratio of ( kg + k, ) to kg for example of 
largest probable error as a function of missile -cone 
half - angle. 
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Cone half- angle, 8 , deg 

Figure 18. - Ratio of (k 3 -k|^) to k 2 for example of 
largest probable importance of kg-k/ as a function of 
missile -cone half -angle. 
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